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Abstract 

I The numerical range of a matrix is studied geometrically via the cone of positive 

semidefinite matrices (or semidefinite cone for short). In particular it is shown that 
the feasible set of a two-dimensional linear matrix inequality (LMI) , an affine section 
of the semidefinite cone, is always dual to the numerical range of a matrix, which is 
therefore an affine projection of the semidefinite cone. Both primal and dual sets can 
also be viewed as convex hulls of explicit algebraic plane curve components. Several 
numerical examples illustrate this interplay between algebra, geometry and semidef- 
' inite programming duality. Finally, these techniques are used to revisit a theorem in 

. statistics on the independence of quadratic forms in a normally distributed vector. 
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1 Notations and definitions 

The numerical range of a matrix A G C"^" is defined as 

W{A) = {w*Aw eC : w eC", w*w = 1}. (1) 

It is a convex closed set of the complex plane which contains the spectrum of A. It 
is also called the field of values, see [SI Chapter 1] and [121 Chapter 1] for elementary 
introductions. Matlab functions for visualizing numerical ranges are freely available from 
[3] and [!!]• 

Let 

A + A* ^ A-A* 
Ao = In, A, = ^—, A2 = ^— (2) 

with In denoting the identity matrix of size n and i denoting the imaginary unit. Define 
J^iA) = {yeFl : F{y) = y^A^ + y^A^ + 2/2^2 >l 0} (3) 
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with >z meaning positive semidefinite (since the Ai are Hermitian matrices, Fiy) has 
real eigenvalues for all y) and denoting the oriented projective real plane (a model 
of the projective plane where the signs of homogeneous coordinates are significant, and 
which allows orientation, ordering and separation tests such as inequalities, see [18] for 
more details). Set J^{A) is a linear section of the cone of positive semidefinite matrices 
(or semidefinite cone for short), see lU Chapter 4]. Inequality F{ii) ^ is called a linear 
matrix inequality (LMI). In the complex plane C, or equivalently, in the affine real plane 
M^, set J^{A) is a convex set including the origin, an affine section of the semidefinite 
cone. 

Let 

p{y) = det{yoAo + yiAi + 7/2^2) 
be a trivariate form of degree n defining the algebraic plane curve 

V = {yeFl : p{y) = 0}. (4) 

Let 

Q = {xeFl : q{x) = 0} (5) 

be the algebraic plane curve dual to V, in the sense that we associate to each point 
y E V a point x E Q of projective coordinates x = {dp{y)/dyo, dp{y)/dyi, dp{y)/dy2)- 
Geometrically, a point in Q corresponds to a tangent at the corresponding point in Q, 
and conversely, see [HI Section V.8] and [TJ Section 1.1] for elementary properties of dual 
curves. 

Let V denote a vector space equipped with inner product (.,.). If x and y are vectors 
then {x,y) = x*y. If X and Y are symmetric matrices, then {X,Y) = trace(X*y). Given 
a set /C in V, its dual set consists of all linear maps from /C to non-negative elements in 
M, namely 

/C* = e V : {x,y)>0, xe /C}. 

Finally, the convex hull of a set /C, denoted conv /C, is the set of all convex combinations 
of elements in )C. 

2 Semidefinite duality 

After identifying C with or P^, the first observation is that numerical range W(74) is 
dual to LMI set J-'{A), and hence it is an affine projection of the semidefinite cone. 

Lemma 1 W{A) = J^{Ay = {{{Aq, W), {Ai, W), {A2, W)) e^l : W e C"^", W h 0}. 

Proof: The dual to J^{A) is 

J'iA)* = {x : (x, y) = (F(|/), W) = j:,{A,, W)yk > 0, ^ 0} 
= {x : Xk = {Ak,W), WhO}, 
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an affine projection of the semidefinite cone. On the other hand, since w*Aw = w*Aiw + 
i w*A2W, the numerical range can be expressed as 

W{A) = {x = {w*Aow, w*Aiw, w*A2w)} 

= {x : Xk = {Ak, W), Who, rank W = 1}, 

the same affine projection as above, acting now on a subset of the semidefinite cone, 
namely the non-convex variety of rank-one positive semidefinite matrices W = ww*. Since 
w*Aow = 1, set W{A) is compact, and conv>V(A) = J'(A)*. The equality W{A) = J^{A)* 
follows from the Toeplitz-Hausdorff theorem establishing convexity of W(A), see fr2{ 
Section 1.3] or [8, Theorem 1.1-2]. □ 

Lemma [T] indicates that the numerical range has the geometry of planar projections of the 
semidefinite cone. In the terminology of [U Chapter 4], the numerical range is semidefinite 
representable. 

3 Convex hulls of algebraic curves 

In this section, we notice that the boundaries of numerical range W(A) and its dual LMI 
set J-" (A) are subsets of algebraic curves V and Q defined respectively in (jl]) and (I5l), and 
explicitly given as locii of determinants of Hermitian pencils. 

3.1 Dual curve 

Lemma 2 J^{A) is the connected component delimited by V around the origin. 

Proof: A ray starting from the origin leaves LMI set J^{A) when the determinant p{y) = 
det J2k Vk^k vanishes. Therefore the boundary of J^{A) is the subset of algebraic curve 
V belonging to the convex connected component containing the origin. □ 

Note that V, by definition, is the locus, or vanishing set of a determinant of a Hermitian 
pencil. Moreover, the pencil is definite at the origin so the corresponding polynomial p{y) 
satisfies a real zero (hyperbolicity) condition. Connected components delimited by such 
determinantal locii are studied in [Q|, where it is shown that they correspond to feasible 
sets of two-dimensional LMIs. A remarkable result of |S] is that every planar LMI set 
can be expressed this way. These LMI sets form a strict subset of planar convex basic 
semi-algebraic sets, called rigidly convex sets (see [9] for examples of convex basic semi- 
algebraic sets which are not rigidly convex). Rigidly convex sets are affine sections of the 
semidefinite cone. 

3.2 Primal curve 

Lemma 3 VV(v4) = conv Q. 



Proof: From the proof of Lemma [Tj a supporting line {x : ^j^XkUk = 0} to VV(v4) has 
coefficients y satisfying p{y) = 0. The boundary of >V(v4) is therefore generated as an 
envelope of the supporting lines. See [16], [HI Theorem 10] and also [HI Theorem l.Sj.D 

Q is called the boundary generating curve of matrix A in [Hj. An interesting feature 
is that, similarly to V, curve Q can be expressed as the locus of a determinant of a 
Hermitian pencil. In the case Q is irreducible (i.e. polynomial q{x) cannot be factored) 
and V is not singular (i.e. there is no point in the complex projective plane such that 
the gradient of p{x) vanishes) then q{x) can be written (up to a multiplicative constant) 
as the determinant of a symmetric pencil, see [Gj Theorem 2.4]. Discrete differentials and 
Bezoutians can also be used to construct symmetric affine determinantal representations, 
see [lOl Section 4.2]. Note however that the constructed pencils are not sign definite. 
Hence the convex hull W(A) is not a rigidly convex LMI set, it cannot be an affine section 
of the semidefinite cone. However, as noticed in Lemma [H it is an affine projection of the 
semidefinite cone. 



4 Examples 



4.1 Rational cubic and quartic 
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Figure 1: Left: LMI set J-" (A) (gray area) delimited by cubic V (black). Right: numerical 
range >V(/1) (gray area, dashed line) convex hull of quartic Q (black solid fine). 
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and 

p{y) = {yo - yi){yo + yif - y^yl 

defines a genus-zero cubic curve V whose connected component containing the origin is 
the LMI set J^{A), see Figure [H With an ehmination technique (resultants or Grobner 
basis with lexicographical ordering), we obtain 

q{x) = Ax\ + 32x2 + ISx^Xg — 1^XqXix\ + Axqx\ — 27x1x1 

defining the dual curve Q, a genus-zero quartic with a cusp, whose convex hull is the 
numerical range W(A), see Figure [H 

4.2 Couple of two nested ovals 




Figure 2: Left: LMI set J-'{A) (gray area) delimited by the inner oval of quartic V (black 
line). Right: numerical range W(A) (gray area) delimited by the outer oval of octic Q 
(black line). 

For 

"0 2 1 + 2i " 
. _ 10 

~ i i 

+ i i 

the quartic V and its dual octic Q both feature two nested ovals, see Figure [2l The inner 
oval delimited by V is rigidly convex, whereas the outer oval delimited by Q is convex, 
but not rigidly convex. 



Figure 3: Left: LMI set J-'{A) (gray area) delimited by the inner oval of quartic V (black 
line). Right: numerical range W(A) (gray area) delimited by the outer oval of twelfth- 
degree Q (black line). 

4.3 Cross and star 

A computer-generated representation of the numerical range as an enveloppe curve can 
be found in [8], Figure 1, p. 139] for 



A = 
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We obtain the quartic 

p{y) = ^(64|/o^ - 52ylyl - b2ylyl + yl + My\yl + yt) 
and the dual twelfth-degree polynomial 

q[x) = 5184x^2 - 299520xJ°a;2 - 299520xJ°x2 + 1954576a;^x^ 

+16356256xgx?x^ + 1954576x^x^ - 5375968x[;x? - 79163552x^x^x1 
-79163552x^x?x^ - 5375968x[;x^ + 7512049x^xf + 152829956x^x?x^ 
-2714586x^xfx| + 152829956x^xfx^ + 7512049x^x^ - 5290740xgxJ° 
-136066372xgxfx^ + 232523512xgxfx^ + 232523512xgxfx^ - 136066372xgx?x^ 
-5290740x2x^° + 1498176xf + 46903680xi°x2 - 129955904x?x| 
-M86148096x?x^ - 129955904xtx^ + 46903680x?x^° + 1498176x^2 

whose corresponding curves and convex hulls are represented in Figure El 
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Figure 4: Left: LMI set J^{A) (gray area) intersection of cubic (black solid line) and conic 
(gray line) LMI sets. Right: numerical range W(A) (gray area, black dashed line) convex 
hull of the union of a quartic curve (black solid line) and conic curve (gray line). 

4.4 Decomposition into irreducible factors 

Consider the example of [H Figure 6, p. 144] with 

1 1 
000001000 
000000100 
000000010 
A= 000000001 
000000000 
000000000 
000000000 
000000000 

The determinant of the trivariate pencil factors as follows 

1 



p{y) 



256 



{Ay^ - 3yoyj - 3yoyl + yf + yiyl){4y^ 



yl 



2\3 



yl) 



which means that the LMI set J^{A) is the intersection of a cubic and conic LMI. 
The dual curve Q is the union of the quartic 

Q 1 = {x : - Sxlxi - ISxlxl - ISxlxl + 27x\ + Mxlxl + 27^2 = 0}, 

a cardioid dual to the cubic factor of piy), and the conic 

g 2 = {x : Xq - 4xi - 4x2 = 0}> 



a circle dual to the quadratic factor of p{y). The numerical range W(A) is the convex 
hull of the union of conv Q i and conv Q 2 , which is here the same as conv Q 1 , see Figure 

m 
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4.5 Polytope 




Figure 5: Left: LMI set ^-"(74) (gray area) intersection of four half-planes. Right: numerical 
range W(A) (gray area) convex hull of four vertices. 

Consider the example of [8], Figure 9, p. 147] with 

" 4 -1 " 
-14 
0-140' 
0-14 

The dual determinant factors into linear terms 

Piy) = ivo + 5yi){yo + Syi){yo + ^yi + ?/2)(yo + 4?/i - 7/2) 

and this generates a polytopic LMI set J^{A) = {y : yo+5yi > 0, yo+3yi > 0, j/o+4|/i+y2 > 
0, yo + 4?/i - 1/2 > 0}, a triangle with vertices (1, -i, 0), (1, -|, |) and (1, -i, -|). The 
dual to curve V is the union of the four points (1,5,0), (1,3,0), (1,4,1) and (1,4,-1) 
and hence the numerical range W(A) is the polytopic convex hull of these four vertices, 
see Figure [51 

5 A problem in statistics 

We have seen with Example 14.51 that the numerical range can be polytopic, and this is 
the case in particular when A is a normal matrix (i.e. satisfying A* A = AA*), see e.g. 
[m Theorem 3] or [SI Theorem 1.4-4]. 

In this section, we study a problem that boils down to studying rectangular numerical 
ranges, i.e. polytopes with edges parallel to the main axes. Craig's theorem is a result 
from statistics on the stochastic independence of two quadratic forms in variates following 
a joint normal distribution, see [1] for an historical account. In its simplest form (called 
the central case) the result can be stated as follows (in the sequel we work in the affine 
plane yo = 1): 
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Theorem 1 Let Ai and A2 he Hermitian matrices of size n. T/ien det(/„+?/iy4i+?/2^2) = 
det(/n + UiAi) det(/n + 2/2^2) if and only if A1A2 = 0. 

Proof: If A1A2 = then obviously det(/„ + yiAi) det(/„ + 2/2^2) = det((/„ + yiAi){In + 
= det(/„ + yiAi + 7/2^2 + 2/iy2^i^2) = det(/n + yiAi + 7/2^2)- Let us prove the 
converse statement. 

Let aifc and a2k respectively denote the eigenvalues of Ai and A2, for k = 1, . . . ,n. Then 
p{y) = det(/„+?/iAi + 2/2^2) = det(/„+|/iy4i) det(/n + 2/2^2) = Ylki'^+Viaik) nfc(l+2/2«2fc) 
factors into linear terms, and we can write p{y) = Ylki^ + + 2/202^) with 01^02^ = 
for all = 1, . . . , ra. Geometrically, this means that the corresponding numerical range 
W(A) for A = Ai + iA2 is a rectangle with vertices (min^ a^, min^ bk), (min^ a^, max^ bk), 
(maxfc Ok, miufc bk) and (max^ a^, max^ bk). 

Following the terminology of [15], Ai and A2 satisfy property L since yiAi + 1/2^2 has 
eigenvalues yiaik + y20-2k for A; = 1, . . . , n. From [151 Theorem 2] it follows that A1A2 = 
A2A1, and hence that the two matrices are simultaneously diagonalisable: there exists a 
unitary matrix U such that U*AiU = diag^^ai^ and U*A2U = diagfca2fc- Since 01^02^ = 
for all k, we have J2k «u«2fc = U*AiUU*A2U = U* A1A2U = and hence A1A2 = 0. D 

6 Conclusion 

The geometry of the numerical range, studied to a large extent by Kippenhahn in [H] - 
see [20] for an English translation with comments and corrections - is revisited here from 
the perspective of semidefinite programming duality. In contrast with previous studies 
of the geometry of the numerical range, based on differential topology [13], it is namely 
noticed that the numerical range is a semidefinite representable set, an affine projection 
of the semidefinite cone, whereas its geometric dual is an LMI set, an affine section of 
the semidefinite cone. The boundaries of both primal and dual sets are components of 
algebraic plane curves explicitly formulated as locii of determinants of Hermitian pencils. 
The geometry of the numerical range is therefore the geometry of (planar sections and 
projections of) the semidefinite cone, and hence every study of this cone is also relevant 
to the study of the numerical range. 

The notion of numerical range can be generalized in various directions, for example in 
spaces of dimension greater than two, where it is non-convex in general [5]. Its convex 
hull is still representable as a projection of the semidefinite cone, and this was used exten- 
sively in the scope of robust control to derive computationally tractable but potentially 
conservative LMI stability conditions for uncertain linear systems, see e.g. [IT]. The 
numerical range of three matrices is mentioned in [I2| Section 1.8]. In this context, it 
would be interesting to derive conditions on three matrices Ai, A2, A^ ensuring that 
det(J„ + yiAi + 1/2^2 + = det(J„ + yiAi) det(4 + Z/2^2) det(4 + Z/s^s)- Another 

extension of the numerical range to matrix polynomials (including matrix Pencils) was 
carried out in [2] , also using algebraic geometric considerations, and it could be interesting 
to study semidefinite representations of convex hulls of these numerical ranges. 

The inverse problem of finding a matrix given its numerical range (as the convex hull 
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of a given algebraic curve) seems to be difficult. In a sense, it is dual to the problem 
of finding a symmetric (or Hermitian) definite linear determinantal representation of a 
trivariate form: given p{y) satisfying a real zero (hyperbolicity) condition, find Hermitian 
matrices such that p{y) = det{^j^ykAk), with Aq positive definite. Explicit formulas 
are described in [9] based on transcendental theta functions and Riemann surface theory, 
and the case of curves {y : p{y) = 0} of genus zero is settled in [10] using Bezoutians. A 
more direct and computationally viable approach in the positive genus case is still missing, 
and one may wonder whether the geometry of the dual object, namely the numerical range 
conv{x : q{x) = 0}, could help in this context. 
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